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Abstract 

Simple toy model is used in order to exhibit the technique of extracting the non-perturbative 
information about Green's functions in Minkowski space. The effective charge and the dynamical 
electron mass are calculated in strong coupling 3+1 QED by solving the coupled Dyson-Schwinger 
equations for electron and photon propagators. The minimal Ball-Chiu vertex was used for simplicity 
and we impose the Landau gauge fixing on QED action. The solution obtained separately in Euclidean 
and Minkowski space were compared, the latter one was extracted with the help of spectral technique. 
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1 Introduction 



In quantum field theory and even in physics at all, the dispersion relations (DRs) were recognized as 
providing useful connections between physical quantities. 

Very recently [1], it was recognized that in the Quantum Chromodynamics (QCD) in Landau 
gauge the full (non-perturbative) gluon propagator is (very likely) an analytic function in the whole 
complex plane of momenta except the positive-timelike [2] real half-axis. This has been achieved by 
the analytical fits of some recent solutions of Dyson-Schwinger equations (DSEs) [3] and by analytical 
parameterization of some contemporary lattice data. The reasonable spacelike domain agreement of 
these fits with the recent Euclidean data give us a good guidance on the possible analytical structure 
of the Green's function in timelike axis of momenta. The similar was argued for the quark propagator, 
however in this case the observed singularities do not occur exactly on - but rather say - very close to 
the real positive half-axis. 

In this paper, we instead to make an analytical guesses of already known Euclidean results we start 
with the assumption of C — Rj^ analyticity of the Green's functions which allows us to solve the DSEs 
system directly in Minkowski space. It is well known that such an assumption is intimately related 
with the formulation of the appropriate dispersion relations (DRs) and spectral representation (SR) 
of Green's functions (for more then the derivation of DRs in the context of QCD see [4]). The main 
advantage of such an approach is the possibility to the knowledge of the propagators at the whole range 
of momentum while the main inconvenience of this method is the necessity of 'inhuman' effort when 
the absorptive parts of Green's functions are actually derived. Unfortunately at the contemporary 
stage of our calculation we are not able to write down the appropriate DRs for ghosts and transverse 
gluons. Due to this fact we confine ourself to the less complicated (at least from the technical point 
of view) case of one flavor strong coupling Quantum Electrodynamics and we leave the involvement 
of Yang-Mills theory for the future study. 

In the present paper and in contrast to [5] wc adopt the Ball-Chiu (BC) vertex [6] which is the 
known minimal ansatz consistent with the Ward-Takahashi identity (WTI). Implementing this into 
the equations for electron and photon wc then solve the corresponding DSE also for the photon 
polarization function. In this sense, the presented study is the extension of the previous numerical 
study of the renormalized electron mass in the strong coupling QED in simplest approximation to the 
DSEs: the bare vertex and bare photon propagators were employed [5]. 

Up to the case of perturbative theory the direct Minkowski space treatment of DSEs is usually 
rather involved and the progress is adequately less when compared with the relatively large amount of 
calculations performed in Euclidean space (for a list of references see the paper [5] and also the paper 
[7] where the approximative analytical solutions are discussed) . Although the spectral approach is 
rapidly getting technically rather involved when one goes beyond lowest order truncation, one of the 
purpose of this paper is to demonstrate that it is still manageable for the considered model with WTI 
respecting vertex. As in the paper [5] we again compare in detail the solutions for fermion propagator 
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in Minkowski and Euclidean space, besides we compare for the first time also the photon propagator. 
There is a another important reason for the use of improved vertex. It is long time known that the 
scalar part of BC vertex dramatically affect the analytical structure of the DSEs [8] (which paper is the 
extended study of the Munczek-Nemirovsky model [9]), for the earlier study of the quark gap equation 
in the axial gauge see [10]. In this place we should stress that the recent study of QCD DSEs confirms 
this suggestion: the bare vertex leads to the complex singularity of the chiral limit quark propagator . 
The authors of Ref. [1] observed that the scalar part of the BC vertex in used plays the crucial role 
in the analytical structure of quark propagator and one can conclude that the inclusion of the scalar 
part of the BC vertex leads again to the real singularity of quark propagator, noting that the same is 
true for quenched QED in the chiral limit (note only that the full Curtis-Pennington vertex [11] was 
used which fact has no large significance due to the Landau gauge employed). 

Contrary to QCD, the considered model here is not an asymptotic free theory and posses the 
additional complication due to the triviality statement. It requires the introduction of ultraviolet 
cut-off function /(A) which well known fact is clearly confirmed by our numerical analysis. However, 
we stress here that QED is trivial and the appropriate solution of DSEs (Euclidean and spectral as 
well as) collapse due to the presence of Landau singularity [12], we regard the strong QED as an useful 
toy model which is clearly reliable when we take cutoff reasonably smaller then the position of the 
expected of Landau singularity. Further, in order to fully specify our model we consider non-zero bare 
electron mass rriQ. 

For comparison we have also calculated the propagators in unquenched approximation (i.e., n(g^) ^ 
0), but with the bare vertex. In that case the vacuum polarization tensor is not gauge invariant and 
contains two independent scalar functions. We impose the transversality by hands before the numerical 
solution. As expected from the earlier studies [13], [14] the results obtained in Landau gauge are very 
close to the results calculated with BC vertex. This pleasant but extraordinary property of Landau 
gauge fixing makes the solutions of DSEs with the bare vertex approximation meaningful. However 
the similar conclusion was made in some QCD studies [15], [16] one should be aware about the inci- 
dence on the analytical structure of the quark propagator. Although when using the bare vertex the 
transverse projection by hands is not theoretically justified and although the analyticity assumptions 
seems not to be fully justified the usage of bare gauge vertices in Landau-like gauges remains popular 
in studies of more complicated gauge models (like extended, walking Technicolor etc. for review see 
[17] ). This is why we believe that showing in detail the effects due to BC vertex is interesting: in 
our studies its effects are less than 10 per cents even for rather large coupling constant. The strong 
coupling QED should serve as an instructive tool for its 'simplicity' and the proposed technique could 
be helpful elsewhere. 

The layout of the article is following: In next section we review the DSE formalism and describe the 
model. The section 3 is devoted to the solution of DESs in Euclidean space. In the section 4 the DSEs 
are written in Minkowski space and the desired DRs for electron selfenergy and gluon polarization 
function are derived. The numerical results are presented in the section 5 and then we summarize. 
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2 Model- Unquenched QED with BC vertex 



In this section we review the basic elements of the considered model and set the notation and conven- 
tions used in the main part of the paper. First of all, let us stress some differences and improvements 
used in this paper when compared with until now published works that are dealing with QED DSEs. 
Using the full photon propagator, instead of the bare propagator that authors of the papers [18], [19], 
[13] (bare and improved vertices, Euclidean space), [5] (bare vertex, Minkowski space) used is an 
important improvement. Further we use the form of the vertex that is consistent with WTI instead 
of bare vertex 7^ that was used in the study of renormalized DSEs in unquenched QED [20]. This 
leads to the solution which reflects not only the effect of running coupling caused by fermion loop 
but also to the vacuum polarization which is automatically transverse. Thus we can call our solution 
as a gauge covariant one since it naturally respect the conservation low that follows from the gauge 
invariance of QED Lagrangian. The three lowest DSEs read: 

[5]-i = [5o]-^-S ; ^ = e'' J ST'^G^.Y, (1) 
[G^u]-' = [GoM.]-' - n^. ; n^u = e^Tr j ST^'SY , (2) 

^^' = Y + eJ sr^'SM = r1 + rt^, (3) 

where Ai is the electron-positron scattering kernel (without annihilation channel) , S is the full fermion 
propagator 

1 F(p2) 

parametrized in a usual way in terms of two scalar functions A,B. Equation (1) then reduces to 
a coupled set of two scalar equations for the functions A, B or equivalently for the dynamical mass 
M = B/A and the electron renormalization function F = A"^. When the interaction is neglected, 
Eq. (4) reduces to the free propagator: Sq^ = i> — 'mo, mo being the bare electron mass. The sixteen 
Lorentz components of Eq. (2) can be reduced to single equation for polarization function 11 



n^-(g) = pA^Vn(g2) ; pA^- = (5M-_:L^), (5) 



by virtue of the gauge invariance q'^H/xi/ = 0. The function Gq" is the quenched approximation (11 = 0) 
to the full photon propagator G>^'^ (q) , which is purely transverse in Landau gauge. 

The functions and F^, in (3) are the longitudinal and transverse parts of the full vertex F'*. 
Multiplying the vertex by the photon momentum p — I one gets the WT identity 

S-\p)-S'\l) = ip-l),T^^{p,l), (7) 

while from its definition Tj,{p, l).{p—l) = 0. Any truncation of DSEs system leading to gauge covariant 
solution of DSEs must involve vertex satisfying WTL Instead of solving own equation for F there 
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exists much economic way. Within the requirement of right Lorentz transformation property, charge 
conservation and the unique hmit when p I (the absence of kinematic singularity) the longitudinal 
part Tj^: 

was found in the paper [6] while the transverse part Tt^{p,l) can be decomposed into the eight com- 
ponent vector basis (see for instance [6] or[21] for the details) 

T!^{p,l) = j2tiTt. (9) 
1=1 

with coefficient functions ti unspecified in general. 

The one loop analysis and determination of ti was performed in the paper [21] in arbitrary covariant 
gauges. The two loop results was obtained in the Feynman gauge also [6]. The minimal version of the 
gauge covariant vertex - BC vertex simply neglect the transverse part and hence is given by the Eq. 
(8). Some better improvement of the full vertex T necessarily differs only by its transverse part r^. 
In chiral symmetric case with given truncation of DSEs and within the requirement of multiplicative 
renormalizability the additional constraints on the transverse pieces of the vertex were found [22]. Also 
the additional information was obtained from the requirement gauge independence of chiral symmetry 
breaking in quenched QED [23]. Although, as it follows form the above two notations the BC vertex 
has not probably most ideal form for some 'definitely conclusive' nonperturbative study of QED, 
nevertheless we use this vertex for its simplicity. The implementation of the methods used in [22], 
[23] and mainly their true Minkowski space extension is far away from triviality and thus remains the 
challenging task for future investigation. 

In the Landau gauge and in the ladder approximation of the electron DSE there is no self-energy 
contribution to the electron renormalization function since the Feynman (F) pole part of photon prop- 
agator Gf^u is exactly canceled by the contribution of longitudinal (LO) Qi^qiy part of G^i,. Explicitly 
their absorptive parts satisfy: 

QAf{u) = -^Alo{uj) = Jda(l- ^ j a^(a)e(a; - a) , (10) 

which leads to A = Ap + Alq = 1 (here av{a) is Dirac coefficient Lehmann function (14), m is the 
physical (pole) electron mass). It is also well known that the property A 1 persists beyond the 
bare vertex approximation and that in angle approximation [24] A = 1 is even valid exactly. The 
above mentioned study [20] of unquenched QED shows rather small and irrelevant violations from 
the identity A = 1 which then turns to few percentage error only when momentum approaches the 
ultraviolet cutoff A (note, the condition ^(0) = 1 was exactly imposed in the paper [20]). Further 
study on dynamical mass generation in unquenched supercritical QED [14], [24] confirm this also in 
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this case. They justify our approximation with at most ten percentage deviation in the infrared region 
(when ^(A) = 1 condition was imposed by the authors. To reduce the complexity of our DSEs we 
explore this nice Landau gauge property and we put explicitly A = I for all momenta. We should 
stress here that this neglection has no effect on the gauge invariance of polarization tensor, i.e. we 
still have 



since the WTI (7) is not violated. 

When we renormalize we adopt the standard notation for the renormalization constants Zi, Z2 
and Z3 (see e.g., [25]). From the approximation employed it follows that Zi = Z2 = 1 which is 
in agreement with the multiplicative renornializability and WTI. Furthermore, the unrenormalized 
vacuum polarization n(^) should be absorbed into the renormalization constant Z3. Similarly, the 
unrenormalized electron self-energy TrS(//)/4 is absorbed into the constant Zjn- 

3 Solution of DSEs in Euclidean space 

The DSEs are often solved in Euclidean space after the Wick rotation A;o k\E is made for each 
momentum. Then the loop integrals should be free of singularities and the Green functions are found 
for positive Euclidean momentum k"^ = A;^ + fcl + + ^4- If there is no additional singularity in 
the complex plane of momenta (that would prohibit the validity of the naive Wick rotation), then 
from the solution for some generic function f{p%) one would get the solution for Minkowski spacelike 
momentum f{p\[)\ j? < 0. The solution for timelike momentum would be in principle obtained by 
the analytical continuation of / to the real axis pj^ > 0. In our case, because of QED triviality in 
four dimensions, the assumption of non-singular behavior holds only with the presence of UV cut-off. 
We presume here, that any numerical attempt to avoid UV cut-off implementation would lead to an 
uncontrolled behavior of the Gell-Mann-Low effective charge. The presence of the ultraviolet cut-off 
is required not only due to the inner consistence (the Wick rotation) but also to ensure the numerical 
stability of our calculation. 

Substituting the Ball-Chiu vertex into the DSEs and employing the projection proposed in [26] (and 
successfully used in the papers [14], [3]) we obtain the following coupled DSEs to be solved numerically: 



(11) 




y + M^y) 



y 




+ 



{M{y)-M{z)) 



y- z 



+ 3M{y){M{y) - M{z)) , 



(12) 



M{x) 
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z{l-Uiz)) 



ZMiy) - iM{y)-M{z))2yxsm^e 
y- z z 



(13) 



where a = c'^/Att and variables x,y,z represent squares of Euclidean momenta, z = x + y — 2-yJyxcosO. 
For details of the derivation of (12), (13) we refer to refs. [26], [14]. 

The gauge invariance q^Ufxu = implies that the photon polarization tensor has to be of the form 

(5. 

In the DSE's formalism the gauge invariance of 11^,^ follows from the gauge covariance of the Ball- 
Chiu vertex. This fact helped the authors of [26] to construct simple recipe how to avoid numerical 
quadratical divergence which would be otherwise presented in equation for IT. Anticipate here that it is 
convenient reduce the photon polarization to a single scalar function also in our Minkowski calculation, 
although there it is not a numerical necessity, but merely matter of technical convenience. 

The first line inside the brackets [. . .] of the Eq. (12) and the first term in the brackets [. . .] 
of Eq. (13) represent the kernels of the bare vertex approximation (with only 7^ retained). They 
give dominant contributions to the dynamical mass of electron and the vacuum polarization as well. 
Neglecting the vacuum polarization effect (putting n(z) = 0) in the Eq. (12), the equation for M 
can be further simplified. This so-called ladder approximation of fermion DSE is represented by one 
dimensional momentum integral equation first derived in ref. [27] and used also in the Euclidean 
confinement study [28]. 



4 Direct treatment in Minkowski space 

Assuming analyticity for complex in C — i?+ Gaussian plain with indicated cut and using the known 
asymptotic behavior of the propagator one can derive the appropriate Lehmann representation (LR) 
for the propagators ( without positivity). The appropriate LR for the fermion propagator in parity 
conserving theory reads 



00 



J p'' — u + ie p — m J p'' — Lo + ie 

where we integrated out the single particle state contribution to the full Lehmann weig 

rS{a — m?). The remaining term cr^c) is assumed to be a real and continuous spectral density that 

originates from the interaction. Similarly we can write for the photon propagator in linear covariant 

gauges 

— h + ie 



where the single photon spectrum rpd{b) can be integrated out as in the previous case. 

Due to the asymptotic the dispersion formula for fermion mass function B requires one subtraction 

. .^(.)^j^a ^'^:'^:^^:!^y (16) 



J q"^ — o + i€ 
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where m(/u) is the renormalized mass at the scale /i. Similar relation could be derived for the function 
A, but here A = 1 and therefore the function B represents renormgroup mass function. 
Renormalized photon polarization function in momentum subtraction scheme reads: 

UnM^J^. ^/^'^'f.-""^ (17) 
J {q^ - uj + ie)[u - n'^) 

where we distinguish two possibly different renormalization scales /x, /tx'. The appropriate renormaliza- 
tion accompanied by the detailed derivation of dispersion relations (16) and (17) is given in the next 
two sections. 

In the perturbation theory the relations for p is usually represented by the series expanded in 
the coupling constant. Here the appropriate relations for p and ps are represented by the integral 
equations involving the Lehmann functions cr's and even the function p itself. Together with two 
additional equations for the a's they form the set of the so called Unitary Equations (UEs) (since 
these are the relations between the imaginary and real parts of propagator (and inverse of propagator) 
functions, here we follow the paper [5] ). In order to derive the UEs recall the well known functional 
identity for distributions 

^ P ■ iirSix' - x) , (18) 



x' — X + ie x' — X 

where P- stands for principal value integration. Making use of the LR for C^^ and of the appropriate 
DR for n in and evaluating the imaginary part of the unit tensor G~^G^"' one arrives to the 
integral equation: 

^7(c)(«) = + K(c) * ' (19) 

where we have adopted a shorthand notation for real functional: 

2 

N ^ , P(a)o-(a;) +a(a)p(x)|^ 

[a*p](a) = P- dx (20) 

Jt a — X 

Notice that the both (19) and (20) are non-zero only for timelike (square of) momentum a (here 
a>T = 4m?). 

For the electron Lehmann weights we get in a similar way: 

f ^ /i + m{p,)f2 m{p)fi + ujf2 

^<'^)^^^= u:-m:\p) ' ^<<^)^^^= u;-mHf,) ' ^^^^ 

/i ^ r + K(c) * Psic)]{^) ; /2 ^ r-^^ + * Ps] (c) , (22) 

which are non-zero only for the timelike u > vr? . 

The physical electron mass is defined by S~^(p = m) = or equivalently M(m) = m. Using the 
dispersion relation (16) the desired relation reads: 
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The residuum value r of pole part of the propagator is fixed already when the renormalization pro- 
cedure is done. Because of our numerical solution is is not convenient to determine r by taking the 
on-shell limit p ^ m directly. The easiest way to evaluate r is the inspection of the real part of the 
identity S~^S = 1 evaluated at some arbitrary scale p^. Choosing for instance p = one gets the 
desired relation 



The original momentum space DSEs are now converted into a coupled set of the UEs (19,21) 
complemented by the subsidiary conditions for the residua and thresholds (23) and (24). 

To solve UEs one has to consider these six integral equations simultaneously at all the positive 
values of spectral variables where the corresponding spectral functions are non-zero. Any internal 
inconsistency (i.e., spacelike Green functions singularities mentioned in the introduction) should be 
seen or felt when the UEs are actually solved. The original momentum space Green functions are then 
obtained through the dispersion relation for the proper function or equivalently by the integration 
of the spectral representation for the full connected propagators S, G. Checking (numerically) this 
equivalence verifies the internal consistence of the method. Compared to the Euclidean approach the 
spectral approach has clear advantage of "already known" analytical continuation at all momenta. 
The disadvantage of the spectral approach is its failure in the (confining) regime where the underlying 
assumptions are not justified. 

4.1 Photon propagator 

In a fixed gauge the photon propagator is fully determined by the gauge independent polarization 
function. We describe below the derivation of once subtracted DR, following from the momentum 
space subtraction procedure for photon polarization tensor. First we briefly review the method in its 
perturbative context. 

In 4 + £ dimensions and for spacelike momentum < the one loop polarization function can be 
written as [29] 



where fifH is t'Hooft dimensionfull scale. The mass-shell subtraction scheme defines Zy, so that 




(24) 



which helps us to avoid dealing with complicated infrared singularities. 




(25) 



jjM/lss^Qj _ Q -^jjig];^ implies that the photon propagator behaves as free one near = 0. Choosing 
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bZj, to cancel entire O(e^) correction we find 

(26) 



127r2 



2 / 
+ 7e — Zn(47r) + In 



and renormalized polarization function in mass-sfiell renormalization prescription satisfies well known 
dispersion relation 

oo 2 

nMA55(^2) ^ n(,2) - iim^n(g2) = / ^^ (^2," + ,,)^ /^('-) (27) 
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with the absorptive part 







7rp(u;) = ^^(l + 2mV'^)\/l-4m2/u;G(tc;- W). (28) 

which is given in many standard textbooks (see for instance [30] , where the result of the integration in 
(27) is written also for timelike momenta). Recall that the one loop 11^"^'^'^ represents also self-energy 
calculated in the popular MS scheme for the special choice of t'Hooft scale Ht'H = "m- [31])- Finally, 
let us remind the definition of the off-shell momentum space subtraction: SZ^ = Ii{iJ,^). Making 
redefinition of the electron charge accompanied by the finite subtraction of (27) we can immediately 
write down the desired dispersion relation (17). 

Now we turn our attention to the derivation of momentum space subtracted Ilji{p^ , /j,^) with the 
dressed propagators and with the full Ball-Chiu vertex included. As mentioned in the previous section 
validity of Ward-Takahashi identity for Fbc naturally leads to the transversality of the polarization 
tensor 

K:Ri<l) = K'q'UuMg') , (29) 
where Vj''^ = {g'^" — ^^) is the transverse projector and capital (R) indicates that renormalized tensor 
(29) must respect gauge symmetry of unrenormalized (U) one. 

The truly massless photons with 11^"^'^'^ (0) = are consequence of the renormalization prescription 
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^u{q') = ^ (30) 

with arbitrary constant C, applied on the full polarization tensor 

U^Jiq) ^ie'J-^ Tr S{1) r^(/, I - q) S{1 - q)] , (31) 
where the explicit dependence of Ball-Chiu vertex on fermionic momenta reads: 

r^(M - q) = 7'- jp^^^imi - qf) - . (32) 

As soon as we use WTI constrained vertex the C independence of resulting 11 is evident but the 
right choice of C facilitates derivation of DR. The reason is that the Pennington-Bloch [26] projector 



T'^Kq) = \ 



(33) 
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cancels the contribution from d+1 space-time metric tensor g^^, which simphfies the actual calculations. 
Let us now derive Hr 

^,22, f 1 ,,,, 

^rKQ ,M ) = |9 ^ m| (34) 

for the case when only part of is retained. Substituting the spectral representation (14) into 
the expression for photon polarization function (34), (31) we immediately get 

3g2^' J (27r)^ 

(tm - 4^) [(/- A)'yv{a) + ^.(g)] [M(fe) + ^.(fe)] 

(/-(?)2-a)(/2-6) • (^^^ 

From now on wc omit the spectral integrals and assume that the presence of any spectral function with 
given arguments automatically implies integration over these variables. Since wc will include explicitly 
the boundaries (thresholds) in step functions in the integral kernels, all integrals can be taken from 
zero to infinity Jq°^. Moreover we label the measure — id^//(27r)^ by di and we also suppress ie factors 
in denominators. Combining the denominators with the help of Feynman parameterization then gives 

n^(,^)(, )-8e JdiJ^ '^\l^ + q^x{l-x)-ax-h{l-x)f ' ^^^^ 
The remaining integral is logarithmic divergent. After the subtraction 

ni?(7,)(9';/^') = ^u{^M) - nc/(7.)(M') (37) 

it leads to the finite dispersion relation. Although this procedure is rather straightforward we present 
for completeness briefly intermediate steps of the derivation. The subtracting procedure (37) yields 
explicitly 



nK(^^)(g2;^2^ ^ gg2| di£dxj^' 



dz 



a,ia)a,{b)[x{l - x)Yiq' - fi')i-2) 



{p + (g2 — ijfi)zx{l — x) + fi'^x{l — x) — ax — 6(1 — x))^ 
- r . Ov{a)ov{h)x{l-x){q'-^P) 

~ (47r)2 io (a;-/x2)(g2_^) ' ^^'') 

^ ' x(l— x) ^ ( / \ -I / 

where (after the loop momentum integration) the substitution z ^ uj = '^i^ made. 

Changing the order of integrations and integrating over x (the appropriate integrals are listed in the 
Appendix A) we obtain the DR: 



e f 

n(7^)(9';/^) = j "'~(^_^2)(^2_^ + ^e) 



AV2(a;,a,6) 



a + b b — a f b — a 
1+ 1 + 



av{a)ay(b) , (39) 
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where A is the weh-known triangle function 



A(x, y, z) = + + — 2xy — 2xz — 2yz . 



(40) 



Considering in the expression above for ay{x) only the delta function parts of spectral functions, i.e., 
rf6{x — m?), we just recover the one loop perturbative result (28) (up to the presence of electron 
propagator residuum r/, which is assumed to be close to 1 when the coupling is small): 



TTPole ( 2 _ 2 «QBd(0) 



Stt 



uj{q^ — UJ + ie) 



1- 



4m^ 



LO 



1 + 



2m^ 



LO 



(41) 



We see immediately that 11^^(0, 0) = as required and that using the projector naturally reproduces 
the perturbation theory in its lowest order. 

Using the prescription (34) we now carry on the derivation for the part of the polarization function 
with the remaining term of Ball-Chiu vertex (second term in rhs of Eq. (32)). First we drop the part 
of the vertex which is proportional to I (because the photon propagator is transverse) and take a trace 
which leads to following finite loop integral: 

-2 >- ^M{l)-M{l-q)' 



ie' 
3^ 



di 



P - {I 



{(7„(a)c7,(6) [2/2 - +31 ■q + 3q^] + ay{b)as{a) [21^ - 8^ - 31 ■ q]} 



(42) 



[{l-qr-aW-b] 

The next intermediate steps of the DR derivation are given in the Appendix B, here we simply present 
the final result: The full polarization function with the Ball-Chiu vertex satisfies the once subtracted 
DR 

(p-tA^) + Prem(t^)) (^^ " f^'^) 



^i^(g^/x2) = I duo 



(q"^ — LO + ie){ijj — ji^ 



(43) 



where follows from (39) and P{rem) from (42) and (69). Explicitly, they read: 



e2 AV2(^,a,6) 

^(7.)(^) = 



LO 



1 + 



a + b b — a 



LO 



LO 



1 + 



b — a 



LO 



ay{a)ay{b)Q (lo - {^/a + Vbf^ 

F{lo, c, c) — F{lo, c, a) + F{lo, b, a) — F(lo, b, c) 



P{rem){^) = ■^(^v{a)as{b)ps{c) 

AV2(a;,c,a) 

t [LO, c, a) = 



(c — 6)(c — a) 



a + c ^ 



00 



LO 



1 + 



uo 



e 



LO 



a + Vcf) 



(44) 



These expressions in their full form have been used in our numerical calculation. No principal value 
integration is necessary and the whole integrand has a regular limit when one spectral variable ap- 
proaches another. Recall that the ordinary integrals over the spectral variables a, 6, c are implicitly 
assumed. The function Trps(c) is simply QM{lo) and it is evaluated in the next section. 
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4.2 Fermion propagator 

In this section we show that the Bah-Chiu vertex 

r^ip -l,p)=r- J^^Z^^Wiip - If) - M(/)] (45) 
substituted to the electron self-energy 

np) = e' j diYS{p-l)Tl{p-l,p)G^,{^ = d,l), (46) 

together with the assumed LR for electron (14) and photon propagator (15) leads to the dispersion 
formula for the dynamical mass (16). In F = 1 approximation we can write 

Tr 

M{p^) = + — S(p) 



= ruo + e — j da . -lY-a ^LiP'^P) p _ i, > (47) 



'Jr y{i>- /)o-i,(q ) + o-s(q)t.m/„ 1 

{p-iy 

where we have adopted conventions and notations of the previous section. Relating bare mass and 
the renormalizcd one by nio = Zm{fi)m{n) and absorbing S(p = n) into the mass renormalization 
constant Zjn{ii) gives the finite mass m(/n) and finite DR for fi independent dynamical mass function 
(16). 

Let us again start with the pure 7^ matrix part of the F^, in S. It leads to the following DR: 



P^ - 1^^ P7m(^) 



-^r-Y i(p ) = / du 

/m2 oj- fi^ p^ ~iu + ie 

A^/\uj, a,b) 



P(>)H = -3(^) a.(a)a,(6)^^^. (48) 

(The derivation is straightforward, see for instance Appendix of Ref. [5]). This is a dominant momen- 
tum dependent part of M. 

Using the remainder terms in Ball-Chiu vertex (45) we get the following contribution to M: 



iv-l? „: 
P 



ps{o)a^{b)a.u{a) 



M^rem){p) = 2e J [(^ _ _ [(^ _ ,)2 _ [^2 _ _ ' ^ 

where we have self-consistently used the formula for difference of the dispersion integrals for M: 

- {p-ir-p^ ^^^^p - ) - M(p )] = y . (50) 

Using the Feynman parameterization (49) is after some algebra transformed into: 

/o "^""l I [{p - ff- L - o{l-lw - oW- hyY ■ ^^^^ 

Matching in (51) two Z-dependent denominators (using a Feynman variables z), making a shift I = l+pz 
and integrating over the momentum I yields the result: 

{Airy Jo Jo Jo Ip'^z^l - z) - axz - o(l - x)z - by{l - z)\ \p^ - o\ 
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which is UV finite by construction. 

It is now easy to write down the DR foUowing from (52) (Some details of its derivation are given in 
the Appendix C). The 'dominant' part following from the pure pole /? = of the photon propagator 
reads explicitly: 



(53) 



TT 2(47r)2 7^2 u — u 

—^ay{a)ps[u)<c) \u — I + — 2— C7^(a)p5(w)e \u — 

To sum it up, the dynamical fermion mass is given by: 

M(p2) = m(/x) + M(^^) (p2) + Mremip") ■ (54) 

Anticipating our numerical results: since the whole Mrem changes the numerical results only slightly 
(as compared to M^^ )), we are approximating its imaginary part in our numerics just by the pole 
contribution MP±. 



5 Numerical solutions and results 

First, let us describe some technical points of our numerical treatment. First, consider the spectral 
approach which is simpler from the numerical point of view. After the formal derivation of the DRs 
for electron self-energy and vacuum polarization function we introduce the positive cut-off A in the 
following way 

F{s) = dx 'P^'l . , (55) 

where the function F represents 11 or S. That is, the absorptive parts of proper function vrp(s) is 
modified by a step function p(s) — >■ p(s)0(A^ — s). The same cut-off is then formally introduced into 
the Lehmann representation for propagators. With the cut-off implemented the set of equations that 
have been numerically solved comprise: the coupled nonlinear integral unitary equations (19), (21) for 
photon and electron Lehmann functions, the equation for absorptive part of self-energies (44,48,74) and 
necessary conditions (23,24). This set of equations is solved by iterations. Then the various Green's 
functions are calculated through the appropriate DRs and considered to be physical for < A^. 

We have found that no other approximation is necessary and the unitary form of DSEs converges 
under the iteration procedure, no matter whether the bare or BC vertex is used. This procedure 
naturally fails when the employed cut-off is rather close to the expected singularity of the running 
constant. Before this numerics fails one observes only a trace of this singularity - the large growth 
of the effective charge at close to A^ Fig. 1). Then the dynamical mass appears to be negative at 
large value of timelike momenta (see the curve for A^ = lO'^ and a(0) = 0.4 in Fig. 4). 

In our Euclidean treatment we adopted the simplest cut-off functions: The Heaviside step function 
Q{—k\ + h?) has been introduced into the kernel of DSEs (12), (13), i.e., upper bound of integrals 
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is replaced oo A|n. The value of the cut-off A^; = 10''M^(0) is taken to be exactly the same as 
in the spectral technique described above, where the zero momentum electron mass M(0) is used 
as a scale. After the subtraction the DSEs for renornialized photon polarization IT and dynamical 
mass M have been solved on suitable grid pf G (0,A^). For the bare vertex the equations can be 
solved by iteration without any numerical problem. The form of the BC vertex makes the numerical 
procedure more difficult, even in the quenched approximation. For the unquenched solution this is even 
more troublesome due to uncontrolled oscillations of the numerical iterations of the running charge 
in the infrared region. Hence we approximate the finite difference in (12), (13) by the appropriate 
differentiation- (the similar trick was used in the paper [3]). Explicitly, we replace: 

(Mfa)-M(.))^ 
y- z 

After making this approximation in the Euclidean DSEs with BC vertex we were able to find the 
numerical solution for the coupling up to the half of the critical coupling for the bare vertex approxi- 
mation. 

The DSEs has to be renormalized with the help of subtraction (both in bare vertex approximation 
and the BC vertex, modified as discussed above). Since in this case the subtraction cannot be done 
analytically (unlike in the Minkowski treatment), one has to implement it numerically, which is not 
straightforward as we describe below. The following expression for polarization function 

nw = /^/(..,)/., (57) 

need to be renormalized so that n/j(0) = (for the explicit form of the function / see (13)). Doing 
this numerically with reasonably high accuracy is not as simple task as for the spectral approach or 
in a perturbation theory. We first solved the equation for the quantity 

fl{x) = jdx [/(x, y) - /(O, y)] /x . (58) 

After finding the solution for 11 we looked the limit x ^ at which 11(0) = K and subtracted this 
constant in order to obtain 'right' 11^^(0;) = Il{x) — K. When one step is not sufficient we repeated 
the procedure. Without this we were not able to find the Euclidean solution an accuracy comparable 
with our Minkowski technique. For instance, for a = 0.2 making the subtraction directly for (57) 
leads to 35 per cent underestimate of a(p^) in the infrared region and 50 per cent underestimate in 
the "ultraviolet" region. Compared to this, five iterations of the procedure described above gives a 
satisfying with 0.2 per cent deviation in the infrared, the subtraction constants in successive inter- 
polations are Ki = -1-0.35, —0.12, -1-0.03, —0.008, -1-0.002. Furthermore, we use the linear interpolation 
and log (perturbative) extrapolation to evaluate functions 11(2;) and M{z) &\, z = x + y — 2xy cos 6 in 
Eqs (12),(13). 

In the both formalisms we use one common renormalization scale (/i^ = 0) to define the running 
coupling 
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where we have exphcitly used a = a(0) and we omit explicit dependence on /i in R-labcl quantities for 
purpose of brevity. We use the same scale in order to renormalize the electron mass. As a mass scale 
of the theory we use M^{p^ = 0) = 1 in arbitrary units. The momentum axis at all figures defined in 
this unit. 

When the coupling a increases the pole mass and the residuum of fermion propagator become 
different from their non-interacting values (r = l;Mp = M(0)). Couple of values of r and Mp 
following from our solutions of DSEs with bare and BC vertex are shown in the Tab.l. The residuum 
r is clearly renormalization (and gauge fixing) scheme dependent quantity. However on the physical 
ground one can expect that 'one particle' contribution to an interacting particle propagator is less 
than one, i.e. we would have naively r < 1, here the values of residua are greater then 1 which is the 
consequence of our subtraction scheme. Not surprisingly, if the similar scheme and the same gauge 
are employed than the property r > 1 survive in the case of quark propagator too [1]. 



a 


0.1 


0.2 


0.4 


Mp/M{0) -BC 


1.044 


1.10 


1.39 


Mp/M(0) -BV 


1.042 


1.09 


1.23 


r -BC 


1.090 


1.22 


1.98 


r -BV 


1.085 


1.19 


1.53 



Tab.l Pole mass and residua of pole part of the electron propagator. The label BC(BV) means the 
results calculated with BC (bare) vertices. The coupling a is the value of running charge at zero 
momenta. 

Let us finally compare our numerical results obtained in the both formalisms with bare or BC 
vertices. The so-called photon renormalization functions G (it is defined by G = a{p'^)/a) are com- 
pared for spacelike momenta in Fig. 1. One sees excellent agreement between the solutions obtained 
in spectral and Euclidean formalism. The one-loop perturbation theory (PT) result is added for com- 
parison. The correct pole mass Mp is used in perturbative formulas. The PT results always below the 
lines corresponding to the DSEs solutions. This can be easily read from the bar vertex form of the 
Euclidean DSE: the decreasing M{x) enhances the function |n|, which being negative must enhance 
a{x). The same functions G are displayed for timelike momenta in Fig. 2., where only the results 
obtained from the unitary equations are presented. Again we would like to stress that the differences 
between the solutions with bare and BC vertices are very small. 

The expected exceptions are solutions with the value coupling constant close to the one for which 
the numerical solutions fail, that is Oc — 0.41 and corresponding G{p^ = — A^) ~ ^G{p^ = A^) ~ 2.6. 
The dynamical fermion mass obtained from DSEs is displayed for spacelike regime of momenta in Fig. 3. 
The small deviation of Euclidean results from the spectral ones can be explained as a numerical error. 
The Fig. 4 shows the difference that follows from the use of different vertices. The absorptive parts of 
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M for bare vertex approximation are also displayed in this Fig. for the same coupling a = 0.1, 0.2, 0.4. 
The negative damp of is observed for a = 0.4. At the end we should remind the reader that the 
bare vertex solution was already obtain by Rakow in the paper [20] . The aforementioned approximative 
independence on the cutoff value was explicitly shown in this work. In order to check the consistence 
we took the cutoff to be as in the paper [20], for instance: M(0)/A = 10~^ together with increasing 
of the renormalized coupling as a = 0.5 then we the agreement between us and the Rakow solution 
was found. 

6 Conclusion 

The Dyson-Schwinger equations for strong coupling QED were solved in the truncation which respects 
the gauge identity. It is the first time when the Minkowski and the Euclidean solutions for lowest QED 
Green functions were made. Working in the Landau gauge we have found a good numerical agreement 
between the solutions obtained in these two technically different frameworks. Prom this we argue 
that the electron as well as the photon propagator posses the standard textbook spectral (Lehmann) 
representations. We showed that at least up to the certain renormalized coupling the proper Green 
functions - the photon polarization function as well as the electron selfenergy- satisfy appropriately 
subtracted dispersion relations. The form of them corresponds with the results already known from 
the conventional perturbation theory. In the other words, there is no significant signal for complex 
singularity of the propagators, which posibility is sometimes sugested in the literature. 

In the other side, the small difference between the Minkowski and the Euclidean solutions carmot 
fully excluded such situation, but from the smallness of observed we can speculate that this effect 
must be rather negligible even for rather strong coupling case a{p^ = 0; A = lO^m) ~ 0.4. 

Triviality of QED was confirmed in our approach. We did not find the possibility to send the 
appropriate ultraviolet cutoff A to infinity with simultaneous keeping the renormalized coupling non- 
zero. With increasing A and with the non-zero bare electron mass we cannot observe a second order 
chiral phase transition as in the paper [20] but instead of this we observe the appearance of Landau 
singularity in the rTinning charge. The obtained results were also compared with bare vertex ap- 
proximation. In that case only the small violation from our 'gauge covariant' solution was observed. 
This difference is still small even for rather strong coupling case a{p^ = 0; A = lO'^m) ~ 0.5. The 
obtained solutions reduce to their perturbative counterpartners in small coupling limit. Furthermore, 
the explicit comparison with the one loop perturbation theory was made. Note, that the triviality 
statement is confirmed in these two later cases too. 

There is still missing link between theory with rriQ = (in chiral symmetry phase as well as in 
chiral symmetry broken phase) and our dispersion technique in used. Nevertheless, we believe that 
the possibility to extract the information about the timelike behavior of Green functions from spectral 
approach remains the main attractive feature when compared with usual Euclidean approach. After 
a certain automatization of the dispersion relations evaluation the method should be extend-able to 
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a more complex theories. The certain progress was already achieved in QCD and the results will be 
published elsewhere. There is also a broad scope for possible future investigations in pure QED:study 
of — matrix property when it is composed from dressed Green- functions, the study of bound states 
with dressed propagators, including transverse correction to the vertex, etc.. 
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A Assorted Integrals 

In this Appendix we list several useful relation. The following integral has been used many times in 
the last step of derivations of DRs: 

Xniuj,a,b)= r e{uj-a/{l-x)-b/x)x''dx , (60) 
Jo 

where a, b are positive real numbers. For the several lowest n it equals: 



Xi{LU,a,b) = n 2 -@{u; - {^/^ + VbY) 



X2{u,a,b) = '-e{oJ-{^/^ + ^/bY), 

/^^/'^{u,a,b)\u! + b- a\\{uj + b- of -2ujh\^, , r- /7^9^ , x 

^3(0;, a, 6) = ^ ' ' ^ ^e(a;-(V^+V6)2), (61) 

where A is the triangle function (40). The variable x in (60) appears from the Feynman parametriza- 
tion of products of the inverse scalar propagators I?i,2: 

D-^D-^ - C dx ^"-Hi-4^-^r(a + /3) 

- Jo [D,x + D2{1 - x)](-+(^)r{a)m ^ ^ 

or from the difference of the propagators I?i,2 

Dr- - or = /' <ix . . (63) 



[Dix + D2{1 - x)](<^+^) ■ 



B Derivation of 11 



R{rem) 



In this Appendix we derive DR for the function Iiu(rem)- To this end we formally interchange the 
labeling of the variables a b in the last term of Eq. (42). Further we substitute I ^ —I + q which 
yields 

2 4e2 f a,(r,)r7,(6)/),s(r) [-1/2 -16^ + 20/ -ry-CK/^] 
^U{rem) ) = ^ j [(/ - g)2 - a] [/2 - 6] [(/ - g)^ - c] [Z^ - c] ' ^^^^ 
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where we have used the dispersion relation formulas for M (16) in order to evaluate their shifted 
argument difference: 

m)-M(i-,)_,^ MO 



P-{l-qY J [{l-qf-cW-c]' 
which is invariant under the shift Z — ^ — Z + g. 

Next it is convenient to rewrite the product of four denominators in (64) making use of 

I{q;a,h) = [[(I - qf - a]f - h]}''. (66) 

It is sufficient to deal only with one term on rhs. of (66), the others are obtained simple by changing 
the spectral variables (the logarithmic divergence appears but it cancels against the same contribution 
of three remaining terms). For instance choosing the variable a, c and making a shift x ^ 1 — x leads 
after the subtraction to: 

4e=^ /•! p a^(a)a,(5)(2 + 4a;-12a;2)(g=^-^2) 

3(47r)2 io 7 -+°^^;;) (a;-/x2)(g2_a;)(c-6)(c-o) " ^ ' 

Integrating over the Feynman variable x yields 

a,(a)a,(5)(2Xo + AX, - UX2){q'' - IJ?) 
3(47r)2 Jo {uj-fj?){q^-u){c-h){c-a) ' ^ ^ 

where Xn is the shorthand notation for the function X„(a;, c, a) introduced in the Appendix A. Gath- 
ering all expressions together one gets the dispersion relation for the polarization function ^R{rem) '■ 

ni^(.em)(9,M) - ^(^J^J-^ (.;-^2)(^2_^) 

J F{u;, c, c) — F{uj, c, a) + F{lo, b, a) — F{u), b, c) 
I ic-b){c-a) 
A^/\u;,c, a) 



F(LO,c,a) = 



2 1 + e(a;- v^+v^ 2'^ (gg 

a; a; V /J ^ ^ 



C Dispersion relation for M^^j-em) 



In this Appendix we derive the absorptive part of Eq. (53). We start from the relation (52) and consider 
the dominant contribution that following from the pole of photon propagator, i.e. cr^{b) = r^S(l3). 
In addition we make the substitution O = "^'^j"!^ which leads to the following double dispersion 
integral: 



(47r)2 7o Jax+o{l-x) Q'^{p'^-Q) - O 

(71) 



2(47r)2 7ii+2 n2(^2_J7) ^2_o 
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The last step is to use algebraic identity 

1 1 r 1 1 



(72) 



{p^ — $7) {p^ — o) — o [p'^ — Q p^ — o ^ 

and re-write the double DR as the difference of single DRs. Substituting (72) into (71) we relabel 
— > a;, o ^ It in the first term and — -u, o — a; in the second one. This cosmetics leads to the 
unsubtracted DR: 

Sp^e^ [ /•°° /■°° {a + u)a^{a)ps{u) „ /■°° , , {u) + a)a^{a)ps{ojY 



P- du 



2(47r)^ [ J2±!i oj-^ioj — u)(j)^ - oj) Jm? Js±i£ u^(u; - u)(p^ - a;) J 

Taking a subtraction at the point p, we get for M(^rem) {P^) o^ice subtracted DR with the weight function 

P{rem) • 

3g2 roo -[ 

= 2(4^^ -X,.^ 

'u + a , , , u^a\ , (j + a , , , , ^ / a; + a\l 
—^a^[a)ps{u)Q I ^ I + —^a^{a)ps{^^)Q \u — I . (74) 
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Figure 1: Charge renormalization function G = a(p^)/a(0) obtained by solutions of DSEs. Each 
beam of hues is labeled by the corresponding coupling a(0). The results of leading order perturbation 
theory (dotted lines) always lie below the DSE result. Only one solution with the BC is shown in the 
figure (for a = 0.4), for smaller coupling the BC solutions would be undistinguishable from the bare 
vertex ones. 
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Figure 2: Spectral solutions for charge renormalization function G for timelike momenta and for 
couplings:a = 0.1,0.2,0.4. The down oriented peaks correspond with the threshold 4M3. 
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Figure 3: The comparison of dynamical mass M{p^) obtained in Euclidean and Minkowski formalism 
for spacelike momenta. For comparison we added also one Euclidean solution calculated with the BC 
vertex (solved with approximations described in the section IV). 
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Figure 4: The absolute values of real and imaginary parts of the electron dynamical mass as obtained 
by solving the unitary equations. The position of the up oriented peaks correspond with the pole mass: 
\Mp, Mp\. The excess for the solution with a = 0.4 is shown, the real part of M becomes negative at 
large momenta. 
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